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Aostract 

An  efficient  parallel  algorithm  to  obtain  maximum 
matcnings  in  convex  bipartite  graphs  i3  obtained.  This 
algorithm  can  De  used  to  obtain  efficient  parallel  algo¬ 
rithms  for  several  scheduling  problems.  Some  examples  are: 
job  scheduling  with  release  times  and  deadlines;  scneduling 
to  minimize  maximum  cost;  and  preemptive  scneduling  to 
minimize  maximum  completion  time. 
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i .  Introduction 


A  convex  bipartite 
l a  , a  ,  ....  a  }  and 
vertices.  £nis  the 
parties : 

(i)  If  (  i ,  j  )  is  an 
and  j-*A; 
two  in  B  • 


graph  G  is  a  triple  ( A,  B,  IS)  .  A  = 
d  =  Id. }  are  disjoint  sets  of 
edge  set.  S^satis?ies  the  following  pro¬ 
edge  of  E,  then  either  i*A  and  j*B  or  i 
i.e.,  no  edge  joins  two  vertices  in  A  or 


(2)  if  (a.,b.)-*E  and  (a.,b.  )*E,  than  (a.  ,  b .  )<E,  l<q<*. 

*1  i  3’’’*  ~ 


Property  (1)  above  is  the  bipartite  property  while  property 
(2)  is  tne  convexity  property.  An  example  convex  bipartite 
graph  is  shown  in  Figure  1.1. 


A 


B 


Figure  A  convex  bipartite  graph. 


F  c  E  is  a  matching  in  the  convex  bipartite  graph 
G=(A,B,£)  iff  no  two  edges  in  F  have  a  common  endpoint. 
Fl={(alb  ),  (a4,b3),  (a^b^Uis  a  matching  in  the  graph  of 

Figure  ?.l  wnile  F2={  ( a,  ,  b^ ) ,  (aj^ ,  b~  )  ,  (a2 ,  ) }  is  not.  F  is 
a  maximum  cardinality  matching  (or  simply  aJ  maximum  match¬ 
ing!  in  G  iff  (a)  F  i3  a  matching  and  (b)  G  contains  no 
matching  tl  such  tnat  |d|>|F|  (|til=number  of  edges  in  H)  . 
The  matching  depicted  by  solid  lines  in  Figure  i.l  is  a  max¬ 
imum  matching  in  that  graph. 

In  what  follows,  we  shall  find  it  convenient  to  have  an 
alternate  representation  of  convex  bipartite  graphs,  it  is 
clear  that  every  convex  bipartite  graph  G»(A,B,E), 
A*{a.  ,  . . .  ,a  j ,  B«*{b.  ,  . . .  ,b  }  is  uniquely  represented  by  tne 
set  5f  triples:  1  m 

T  *  t(i,s. ,h. ) ll<i<n} 
si*min{ j | tai7bj )«  1} 
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h^=max{ j I (ai#  bj ) 4E} 

intuitively,  in  the  triple  representation,  we  expli¬ 
citly  record  the  smallest  (s^)  and  highest  (h^)  index  ver¬ 
tices  to  which  each  a.  is  connected.  For  the*  example  of 
Figure  1.1,  we  have  T*=  {(1,1,2),  (2,3,3),  (3,d,3),  (4,3,3), 
(5,1,3)  j  . 

As  an  example  of  the  use  of  matchings  in  convex  bipar¬ 
tite  graphs,  consider  tne  problem  of  scheduling  n  unit  time 
jobs  to  minimize  the  number  of  tardy  jobs.  in  this  problem, 
we  are  given  a  set,  of  n  jobs.  Job  i  has  a  release  time  r. 
and  a  due  time  d..  It  requires  one  unit  of  processing.  A 
subset  F  of  J  is  a  f easiole  subset  iff  every  job  in  F  can  be 
scheduled  on  one  machine  in  such  a  way  that  no  job  is 
scneduled  before  its  release  time  or  after  its  due  time.  A 
feasible  subset  F  is  a  maximum  feasible  suoset  iff  tnere  is 
no  feasible  subset  U  of  J  such  that  IqI > I jf . 

A  maximum  feasible  subset  F  can  be  found  by  transform¬ 
ing  the  problem  into  a  maximum  matching  problem  on  a  convex 
bipartite  graph.  Without  loss  of  generality,  we  may  assume 
that  all  the  r.s  and  d.s  are  natural  numbers;  min[r. }=J; 
r.<d. ,  l£i£n;  and*max  {d.jin.  The  convex  bipartite  graph 
corresponding  to  J  is  ^iven  by  the  triple  set  T  =  {(i,  s., 
hi)  I  s^=r^, h^=d^-l ) .  Figure  1.2  shows  an  example  joD  set 
and  thexcorresponding  convex  bipartite  graph  G.  Vertex  i  of 
A  represents  job  i  while  vertex  i  in  B  simply  represents  the 
time  slot  [i,i+lj.  There  is  an  edge  from  job  i  to  time  slot 
C j, j+1]  iff  r.£j<d. .  Hence,  every  matching  in  G  represents 
a  feasible  subsetof  J.  Also,  corresponding  to  every  feasi¬ 
ble  subset  of  J  there  is  a  matching  in  G.  Clearly,  a  max¬ 
imum  cardinality  feasible  subset  of  J  can  be  easily  obtained 
from  a  maximum  matching  of  G.  In  addition,  a  maximum  match¬ 
ing  also  provides  the  time  slots  in  which  the  jobs  snould  be 
scheduled. 

We  shall  sea  several  other  examples  of  tne  application 
of  matchings  in  convex  bipartite  graphs  in  later  sections. 

Glover  [4J  has  obtained  a  rather  simple  algorithm  to 
find  a  maximum  matching  in  a  convex  bipartite  graph 
G=*  ( A ,  B ,  i2 ) .  Let  h.=max{jl  (a.,b.)«  E),  l_<^i^  I  A.  I  .  Glover's 
algorithm  considers  the  vertices  in  B  one  by  one  starting  at 
b^ .  We  first  determine  the  set  R  of  remaining  vertices  in  A 
to  which  the  vertex  b.  currently  being  considered  is  con¬ 
nected.  Let  q  be  such  that  a_-*R  and  h  »  min  {hi.  Vertex 

^  ^  a  4R  P 

b.  is  matched  to  a  and  a  deleted  from  th8  grapn.  The  next 
vertex  in  B  is  now  considered.  Figure  1.3  specifies 
Glover's  algorithm  more  formally.  By  using  appropriate  data 
structures,  Glover's  algorithm  can  be  implemented  to  run  in 
0(mn)  time. 
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Figure  1^.2 

procedure  MATCH ( 3 ,  m ,  n ,  h  ) 

// /G-(a,  B,  E.)  is  a  convex  bipartite  graph .  | A I =n ,  I  B I  =m/  / 
//n  =max{ j I (a ., b . )<E} . A  maximum  matching// 

/ / Mis  found//  3 

set  M;  inteqer  m,n,h.  graph  G 

M  <—  <t> 

for  j.  1  to  m  do 

let  R  be  the  set  of  remaining  vertices  to  wnich  b 
is  connected 


if  R/fJ  then  q  «—  i  such  that  j«R  and  min  (h 

J  P 


} 


M  «-  M  U( (a  , b- ) } 
delete  a  from  G 
endif  ^ 


p*R 


end  for 
return  ( ) 
end  MATCH 


Figure  1_.3_  Glover's  maximum  matching  algorithm. 


In  this  paper,  we  are  primarily  concerned  with  parallel 
algorithms.  In  Section  2,  we  obtain  a  parallel  algorithm 
for  maximum  matchings  in  convex  bipartite  graphs.  In  Sec¬ 
tion  3,  we  use  this  parallel  algorithm  to  obtain  efficient 
parallel  algorithms  for  three  scheduling  problems. 

Tne  parallel  computer  model  used  is  tne  shared  memory 
model  (SMM) .  This  is  an  example  of  a  single  instruction 
stream  multiple  data  stream  (SIMO)  computer.  In  a  3 MM  com¬ 
puter,  tnere  are  p  processing  elements  (PEs) .  Each  PE  is 
capable  of  performing  the  standard  arithmetic  and  logical 
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operations.  The  PEs  are  indexed  J,l,...,p-1  and  an  indivi¬ 
dual  PH  may  oe  referenced  as  in  PE(i).  Eacn  PS  Knows  its 
index  and  nas  some  local  memory.  in  addition,  there  is  a 
global  memory  to  which  every  PE  has  access.  The  PEs  are  syn¬ 
chronized  and  operate  under  the  control  of  a  single  instruc¬ 
tion  stream.  An  enable/disable  mask  may  be  used  to  select  a 
subset  of  the  PEs  that  are  to  perform  an  instruction,  only 
the  enabled  PEs  will  perform  tne  instruction.  Disabled  PEs 
remain  idle.  All  enabled  PEs  execute  the  same  instruction. 
The  set  of  enabled  PEs  may  change  from  instruction  to 
instruction . 

If  two  PEs  attempt  to  simultaneously  read  the  same  word 
of  tne  shared  memory,  a  read  conflict  occurs.  if  two  PEs 
attempt  to  simultaneously  write  into  the  same  word  of  the 
shared  memory,  a  write  conflict  occurs.  Throughout  this 
paper,  we  shall  assume  tnat  read  and  write  conflicts  are 
prohibited . 

The  reader  is  referred  to  l2]  for  a  list  of  references 
dealing  with  graph  algorithms,  matrix  algorithms,  sorting, 
scheduling,  etc.  on  a  EMM  computer. 


_2.  Parallel  Matching  In  Convex  Bipartite  Graphs 

In  Section  1,  we  showed  that  every  instance  of  the  problem 
of  scheduling  jobs  to  minimize  the  number  of  tardy  jobs 
could  be  transformed  into  an  equivalent  instance  of  the  max¬ 
imum  matching  in  a  convex  bipartite  graph  problem.  It 
should  be  evident  that  the  reverse  is  also  true.  Hence,  the 
two  problems  are  isomorphic.  A  parallel  algorithm  for  a 
special  case  of  the  job  scheduling  formulation  was  ODtained 
oy  us  in  L1J.  In  this  special  case,  it  was  assumed  that  all 
jobs  nave  tne  same  release  time.  This  corresponds  to  the 
case  when  the  convex  bipartite  graphs  are  of  the  form 
T={ ( i , s^, h^ ) | 1 <i<n}  s^=c,l^i^n  for  some  c. 

We  shall  now  proceed  to  show  how  the  solution  for  the 
special  case  described  above  can  be  used  to  solve  the  gen¬ 
eral  case  when  all  the  r.s  are  not  the  same.  This  will  be 
done  using  the  binary1  tree  method  described  by  Dekel  and 
Sahni  [2J.  Rather  than  specify  the  new  algorithm  formally, 
we  shall  describe  how  it  works  by  means  of  an  example. 

A  convex  bipartite  graph  is  shown  in  Figure  2.1.  For 
this  graph,  |A|*14  and  }al=13.  The  s-  and  h.  values  associ¬ 
ated  with  each  vertex  of  A  are  given  In  the  first  two 
columns  of  this  figure.  The  first  step  in  our  parallel 
algorithm  for  maximum  matching  is  to  sort  the  vertices  in  A 
in  nondecreasing  order  of  s. .  Vertices  with  tne  same  s^  are 
sorted  into  nondecreasing  order  of  h^ .  For  our  example,  the 


B 


Figure  2^.1 

:  of  this  reordering  is  shown  in  Figure  2.2 . 
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Following  the  sort,  we  identify  the  distinct  s.  values. 
Let  tnese  Da  R  ,  R_  ,  .  .  . ,  R  .  Assume  that  R.  <  R .  <  .  .  .  < R.  .  Let 
R,  =max{d.}+l.  For  our  example,  k=4  and^Rt 1 : k+1 )  *  (1,  4, 
9,12 ,  14)* 

We  are  now  ready  to  use  the  binary  tree  method  of  [2]. 
The  underlying  computation  tree  we  shall  use  is  the  unique 
complete  binary  tree  with  k  leaf  nodes.  Figure  2.1  shows 
tne  complete  binary  trees  with  4,  5,  and  6  leaf  nodes.  For 
our  example,  k=4  and  we  use  the  tree  of  Figure  2.3(a).  With 
each  node,  P,  in  the  computation  tree,  we  associate  a  con¬ 
tiguous  subset  {u, u+1, u+2, . . . , v}  of  the  vertices  in  B.  This 
subset  is  denoted  [u,v].P  or  simply  [u,v]. 


(a)  4  leaves 


Figure  2.3^t  Complete  binary  trees. 


Let  the  leaf  nodes  of  the  computation  tree  be  numoered 
1  through  k,  left  to  right.  If  P  is  the  ith  leaf  node,  then 
Cu,v].P  is  LRj* Rj. i J  (i.e.,  u=R.  and  v=R.+.-l).  if  P  is 
not  a  leaf  riodeTthen  the  subset  of  B  associated  with  P  is 
Lu,vJ.LC(P)  U  [u,v].RC(P)  where  LC(P)  and  RC(P)  are,  respec¬ 
tively,  the  left  and  right  children  of  P.  The  subsets  of  B 
associated  with  each  of  the  vertices  in  the  computation  tree 
for  our  example  are  shown  in  Figure  2.4.  The  number  in  each 
node  of  this  tree  i3  its  index. 


Let  P  be  any  vertex  of  the  computation  tree.  Let  [u,vj 
be  the  subset  of  B  associated  with  P.  The  subset  of  A  avail¬ 
able  for  matching  at  node  P  is  denoted  M(P)  and  is  defined 
to  be: 


iv  & 

[4,8]  [9,11] 


[12,13] 


Figure  2.4 


M(P) =[i | u£s^<v ] 


For , example, 

M(l)=[l,2,  .  ...14); 
M(2)=[l,2,3,4,5,6,7,12); 
M(4)={2, 3,4,7}; 
etc . 


The  subset  M(P)  of  A  vertices  available  for  matching  at 
P  may  be  partitioned  into  three  subsets  Q(P),  I(P),  and 
T ( P ) .  0(P)  is  a  maximum  cardinality  subset  of  M(P)  that  may 
be  matched  with  vertices  in  [u,v].P  by  algorithm  MATCH;  tnis 
subset  is  called  the  matched  set .  I(P)  denotes  the  infeasi¬ 
ble  set .  It  consists  of  all  vertices  i-*M(P)-Q(P)  such  that 
h-<y.  The  transferred  set  T(P)  consists  of  all  vertices 
iiM(P)— d(P)  such  that  h^>v. 

Consider  node  2  of  Figure  2.4.  The  matching  problem 
defined  at  this  node  is  given  in  Figure  2.5.  If  Glover's 
algorithm  is  used  on  this  graph,  then  [1,2,3,4,5,7,12] 
defines  a  subset  of  A*  that  can  be  matched  with  vertices  in 
B 1 .  Further,  this  gives  a  maximum  matching.  Hence, 
Q(2)={1,2,3,4,5,7,12};  1(2)  =6;  and  T(2)=(6.  Observe  that 
|Q(1 ) I  is  tne  size  of  a  maximum  matching  in  the  original 
convex  bipartite  graph.  Also,  T ( 1 )  and  I(1)=A-Q(1). 

We  shall  make  two  passes  over  the  computation  tree. 
The  first  pass  begins  at  the  leaves  and  moves  towards  the 
root.  During  this  pass,  the  Q,  I,  and  T  sets  for  each  node 
are  computed.  The  second  pass  starts  at  the  root  and 
progresses  towards  the  leaves.  In  this  pass,  the  Q  set  for 
each  node  is  updat  @d  so  as  to  correspond  to  the  set  of  A 
vertices  ma  -!hed  by  f  .over's  algorithm  to  the  B  vertices 
associated  ' tb  node. 
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Pass  _I 

in  this  pass,  we  make  extensive  use  of  the  parallel  algo- 
ritnm  developed  in  CU  for  the  case  when  all  the  s.s  are  the 
same.  For  our  purposes  here,  it  is  sufficient  to  Know  the 
sequential  algorithm  (FEAS  of  C 1 3 )  that  this  parallel  algo¬ 
rithm  is  based  on.  This  sequential  algorithm  is  given  in 
Figure  2.6.  For  convenience,  this  has  been  translated  into 
the  graph  language  used  here. 

line  procedure  FEAS(n,u,v) 

// find  a  maximum  matching  of  vertices  in  A  onto  the  B 
3et  [u,v]  for  every  vertex  i*A,  s.*u // 

1  global  MAT ( 1 : n ) ;  set  A;  integer  n,i,v,i,j 

2  sort  A  into  nondecrea3ing  order  of  h. 

2  MAT(lsn)  «—  J/ /initialize//  1 

4  j  u 

5  for  i  < —  1  to  n  do 

6  case 

7  :j>v:  return(j)  //all  vertices  in  B  matched// 

8  s j£h . s  //select  i//  j  «—  j+1,  MAT(i)  1 

9  end  case 
l<d  end  for 

11  return ( j) 

12  end  FEAS 

Figure  2 . 6 


An  axaminaaon  of  Glover's  algorithm  (Figure  1.3) 
reveals  that  it  performs  exactly  as  does  procedure  FEAS  when 
the  restrictions  and  simplifications  applicable  to  FEAS  are 
incorporated  into  Figure  1.3. 

Hence,  for  a  leaf  node  of  the  computation  tree,  the  Q 
set  may  be  obtained  by  a  direct  application  of  procedure 
FEAS  (or  its  parallel  equivalent) .  For  example,  for  node  4 
of  Figure  2.4,  we  nave  A  =  M(4)  =  {2, 3, 4, 7};  h2=4?  h  =3; 
n  =4;  'n  =2;  u=l;  and  v=4.  Using  FEAS,  we  ootain 
dT4 )= { 7 , i , 2 } .  dote  tnat  J  consists  of  exactly  those  ver¬ 
tices  i  with  MAT(i)=l.  I()  consists  of  exactly  those  ver¬ 
tices  i  witn  MAT (  i )  =(6  and  h._<v.  The  remaining  vertices  form 
T ( ) •  The  matched  set  Q,  transferred  sec  T,  and  infeasible 
set  I  for  eacn  of  the  leaves  in  our  example  are  shown  m 
Figure  2.6.  Hull  sets  are  not  explicitly  snown.  So,  for 
node  4,  i(4)=**6;  T(4)={4}?  and  0(4  )  =  (  7 , 3 , 2  }  .  The  sets  are 
ordered  by  h^. 

For  a  non leaf  node  P,  tne  Q,  1,  and  T  sets  may  be 
obtained  by  using  the  Q,  I,  and  T  sets  of  tne  children  of  P. 
Let  L  and  R,  respectively,  be  the  left  and  right  children  of 
p.  To  determine  Q(P),  we  use  procedure  FEAS  with  u=uR  and 
v=v  (iu  ,  v  j  is  the  subset  of  B  associated  witn  the  right 
child  RH  or  P) .  The  A  sat  consists  of  T(L)  U  0(R).  Since 
both  T(L)  and  Q(R)  are  already  sorted  by  h^,  the  sort  of 
line  2  of  FEAS  can  be  replaced  by  a  merge. 

Let  S  be  the  subset  of  T(L)  U  Q(R)  that  has  MAT()=1 
following  the  execution  of  FEAS.  The  following  theorem 
establishes  that  Q(L)  U  S  is  a  maximum  cardinality  subset  of 
M(P)  that  may  be  matched  with  vertices  in  [n,vj.p.  Hence, 
G(P)=Q(L)  U  S. 

Tneorem  2.1_:  S  as  defined  above  is  a  maximum  cardinality 
subset  of  M(P)  that  may  be  matcned  with  vertices  in  [u,v].P 
using  algoritnm  MATCH. 

Proof :  The  proof  is  by  induction  on  the  distance  of  P  from 
the  farthest  leaves  in  the  subtree  of  wnicn  it  is  a  root, 
if  this  distance  is  1,  then  Q(L)  and  Q(R)  are  maximum  cardi¬ 
nality  subsets  of  M(L)  and  M(R)  tnat  can,  respectively,  be 
matched  by  Glover's  algorithm  with  vertices  in  [u,vJ.L  and 
Cu,v].R.  If  tnis  distance  is  greater  than  1,  then  0(b)  and 
U(R)  satisfy  this  maximum  cardinality  matching  requirement 
by  induction. 

As  far  as  node  P  is  itself  concerned,  we  see  that  only 
vertices  in  M(L)  are  candidates  for  matching  with  vertices 
in  Cu, ,v^J  (recall  that  for  vertices  in  M(R),  the  value 
exceeds  v  ) .  Furthermore,  when  Glover's  algorithm  is  used 
with  the  A  set  being  M(P)  and  the  B  set  being  = 
£u,vj.P,  vertices  in  B  are  considered  in  the  order  u^, 
u.+l, ...,vr ,  u_,...,vD.  Hence,  Q(L)  is  precisely  the  subset 

Li  Li  K 


il 


Figure  2.7^:  Results  of  first  pass, 
of  M(P)  chat  gets  matched  with  vertices  in  [u^.v^]. 


The  candidates  for  the  remaining  vertices  in  B,  i.e., 
£u  ,vRJ  are  clearly  T(L)  U  M(R).  From  the  way  Glover's 
algoritnm  works,  it  is  also  clear  that  the  vertices  of  T(L) 
U  M(R)  that  will  get  matched  to  Cu_,vR]  are  a  subset  of  f(L) 
J  Q(R).  Let  this  subset  be  S' .  We  wish  to  show  that  S  is  a 


legitimate  choice  for  S'.  First,  we  show  that  S  represents 
a  feasible  matching.  Then,  we  shall  show  that  S  is  in  face 
selectable  by  Glover's  algorithm. 


V»e  know  that  Q(R)  can  be  matched  into  £uR,vR].  Let  ^ 
any  such  matching.  Since  S  is  selected  by  FEAS,  we  know 
that  every  vertex  in  S  can  be  matched  to  a  vertex  in  [uR,vR] 
in  a  3uch  a  way  that  no  vertex  j  in  S  is  matched  to  a  vertex 
with  index  greater  than  h ^ •  Consider  a  matching  W  tnat 
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meets  this  condition.  Now  suppose  that  soma  vertex  j  in  S 
is  matcned  to  a  vertex  q  in  LuRlv„]  with  index  less  than  s.. 
Clearly,  j  mu3t  be  a  member  ofKQ(R)  (as  all  vertices  in  T(£) 
nave  an  s  value  less  than  uR) .  Suppose  that  j  is  matched  to 
in  Z.  If  j^,  is  free  in  W,  then  the  matching  of  j  in  W 
may  be  changed  from  q  to  j. .  If  j.  is  not  free,  then  sup¬ 
pose  it  is  matched  to  j_.  If  j2*T(L),  then  j_  may  be 
matched  to  q  and  j  to  j.  (note  that  q<g.  and  so  s.  *q<h  .  ). 

if  j24Q(R)»  then  suppose  that  is  matcned  to  in  Z.  It 
is  easy  to  see  that  If  q=j^  or  is  free  in  W,  then 

we  may  matcn  j  to  ana  j2  to  q.  Otherwise,  let  be  the 
vertex  matched  to  in  W .  it  should  be  clear  that  we  can 
continue  in  this  way  and  modify  W  so  tnat  j  is  matched  to 
jl»  to  j^,  j.  to  jc»  etc.  In  tne  new  matching,  there  is 
one  fewer  vertex  or  S  that  is  matched  to  a  vertex  with 
smaller  s  value. 

Repeating  tne  aoove  construction  several  times,  W  can 
be  transformed  into  a  matching  such  that  every  vertex  j*S  is 
matched  to  a  vertex  q  such  that  uR£s  £q^.h  .£vR.  Hence,  5 
represents  a  feasible  matching.  ^  ^ 

Let  S'  (as  defined  earlier)  be  the  subset  of  Q(R)  U 
T(L)  matched  by  Glover's  algorithm  to  the  vertex  3et 
[uR,vR].  We  shall  now  proceed  to  show  that  5  is  a  valid 
choice  for  S'.  Let  Z  be  any  matcning  of  Q(R)  into  lu  ,v  ]. 
Let  Y  be  a  matching  of  S*  in  which  all  vertices  in  Q(RT  nKS’ 
are  matched  to  the  same  vertex  in  [uR,vR]  as  in  tne  matching 
Z.  Let  W  be  a  corresponding  matching  for  S.  The  existence 
of  the  matchings  Y  and  W  is  a  consequence  of  the  construc¬ 
tion  used  to  show  the  feasibility  of  S. 

From  the  definition  of  S’,  it  follows  that  5‘^S.  Also, 
from  the  working  of  FEAS,  it  follows  that  S#5S*  .  Let  j*S  be 
a  vertex  with  least  h .  such  that  j  /  S'  .  If  no  sucn  j 

exists,  then  S=S 1 .  Assume  that  j  is  matched  to  q  in  W.  If 

q  is  free  in  Y,  then  S'  cannot  be  of  maximum  cardinality. 
So,  let  P*S '  be  matched  to  q  in  Y.  By  definition  of  Y  and 
W,  p  /  S.  Also,  from  the  order  in  which  FEAS  considers  ver¬ 
tices,  h  >n..  Hence,  S'u(j)-{p}  is  also  a  subset  selectaole 
by  GlovarTs^algorithm.  S'  U  (j}-tp)  agrees  with  3  in  one 

place  more  than  does  S'. 

By  repeating  this  interchange  process.  S’  may  be 

transformed  into  S  with  the  result  that  S  is  also  a  maximum 
cardinality  subset  of  Q(R)  U  T(L)  that  is  selectable  by 
Glover's  algorithm  for  matching  in  £uR,vR]. 

Hence,  Q(L)  U  S  is  a  maximum  cardinality  subset  of  M(P) 
selectable  by  Glover's  algorithm  for  matching  in  Lu,vJ.P. 
[  ] 


Once  d(P)  is  Known,  T(P)  and  I(P)  are  easily  computed. 
Actually,  as  I(P)  is  never  used,  we  may  omit  its  computa¬ 
tion.  Figure  2.6  snows  tne  0,  X,  and  T  sets  (except  when 
empty)  for  all  nodes  in  our  example. 


Pas3  2 

Xn  the  second  pass,  for  each  vertex  P  of  the  computation 
tree,  we  compute  a  set  a' (P)  which  represents  the  set  of  A 
vertices  matched  toy  Glover's  algorithm  with  the  set  [u,v].P. 
With  repsect  to  the  matching  shown  by  solid  lines  in  Figure 
2.1,  we  see  that  if  P  is  the  root,  then  Q,  (P)  =  [1,  2,  i, 
4,  5,  7,  3,  9,  Id,  11,  12,  14};  if  P  is  node  3  of  the  compu¬ 
tation  tree,  then  Q' (P)  =  [3,  9,  10,  11,  14}. 

Xf  P  is  the  root  node,  then  Q'(P)*Q(P),  by  definition 
of  Q(P) .  Let  P  be  any  nonleaf  node  for  which  Q'(P)  has  been 
computed.  Let  L  and  R  be  the  left  and  right  children. 


that  both  V  and  Q(L)  can  be  maintained  so  tnat  they  are  in 
nondecreasing  order  of  h.  and  that  W  is  also  in  nondecreas- 
ing  order  of  V).  Q' (L)  consists  of  the  first  min  {|W|, 
v.-u.+l}  vertices  in  W.  The  correctness  of  tnis  statement 
may  fee  established  by  induction  on  the  distance  of  P  from 
the  root,  d' (R)  is  readily  seen  to  be  a'(P)-J'(L).  Figure 
2.7,  shows  the  Q ' (  )  sets  for  all  the  vertices  in  the  compu¬ 
tation  tree  of  our  example. 

From  tne  Q' (  )  sets  of  the  leaves,  the  matching  is 
easily  obtained.  Xf  P  is  a  leaf,  and  [u, v J . P=[a, b] ,  then 
tne  first  vertex  in  Q' (P)  is  matched  with  a,  the  second  witn 
a+X,  etc.  (note  that  Q* (P)  is  in  nondecreasing  order  of  h^). 
The  matching  for  our  example  is  also  given  in  Figure  2.7. 


respectively, 
[u  ,vRJ.  Let 
ordered  set 


of  P.  Let  [u,vj 
V  =  { j I j*Q 1 (P)  and 
obtained  by  merging 


L  =  C u , , v  j  and  [u,v].R  = 
h.<u  }V  ^Let  W  be  tne 
together  V  and  0(L)  (note 


Complexity  Analysis 

Tne  initial  ordering  of  A  by  h.  and  within  h.  oy  S.  can  be 
done  in  O(log^n)  time  using  n/2  P£s([6j  and  [7];.  During 
the  first  pas3,  the  computation  of  Q( )  requires  the  use  of 
FEAS  and  a  merge.  The  use  of  FEAS  (without  the  sort)  taxes 
O(logn)  time  and  requires  0( |M(P) |/log|M(P) | )  PEs.  The 
merge  at  node  P  takes  O(logn)  time  with  |M(P)l/2  PEs. 
Since,  Q  can  be  computed  in  parallel  for  all  nodes  on  the 
same  level  of  the  computation  tree,  O(logn)  time  is  needed 
per  level.  The  total  time  for  the  first  pass  is  0(log  n) 
and  n/2  PEs  are  needed.  Pass  2  requires  only  some,merging 
per  node.  The  total  cost  of  this  pass  is  also  O(log^n)  and 
n/2  PEs  suffice. 
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matching={ (7,1) , (3,2) , (2,3) , (4,4) , (1,5) , (5,6) , (12,7) , (9,9) , (8,10) 
(11, 11), (14, 12), (10,13)} 


Figure  2^3  Second  pas3 

Hence,  tne  overall  complexity  of  our  parallel  algorithm 
for  maximum  matching  in  convex  bipartite  graphs  is  0(log  n). 
The  PE  requirement  is  0(n) . 

Another  complexity  measure  often  computed  for  parallel 
algorithms  is  the  effectiveness  of  processor  utilization 
(EPU)  (see  LlJi  C2J*  and  [9]).  For  any  problem  P  and  paral¬ 
lel  algorithm  A,  this  is  defined  as  follows: 

£PU(P,A)  - 

complexity  of  fastest  sequential  algorithm  for  P 
complexity  of  A  *  number  of  PEs  used  by  A 


For  2our  algorithm,,  we  have  an 
O(run/(log^n*n) )  ■  0(m/log^n)  (recall  that  i 


EPU 
m*  |  B I  )  . 


that 
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j_.  Applications  to  Scheduling 

We  nave  already  seen  one  application  of  maccnings  in  convex 
bipartite  graphs  to  scheduling.  This  application  «ras  to  a 
proolem  that  is  actually  isomorphic  to  tne  matching  proDlem; 
viz.  Find  the  maximum  number  of  unit  time  joos  from 
J={(r^,d.)l  l£i£n}  that  can  be  scheduled  on  a  single  machine 
such  th&t  no  job  is  scneduled  before  its  release  time  r.  or 
after  its  due  time  d. .  In  tnis  section,  we  snail  loo*  at 
three  other  scheduling  problems  that  can  be  solved  in  an 
efficient  manner  using  the  parallel  matching  algorithm 
developed  in  previous  section. 


iL'A  Scheduling  to  Minimize  the  Maximum  Cost 


In  this  problem,  we  are  given  n  unit  time  jobs.  Associated 
with  each  job  i3  a  release  time  r.  and  a  nondecreasing  cost 
function  c-(),  l£i<n.  Let  S=(  s,  ,  s^ »  •  •  •  #  3n)  be  a  one  machine 
schedule  for  the  n  jobs .  s .  denotes  the  start  time  of  job 
i.  Let  q.  be  as  defined  below: 


ci<siJ  3iiri 


oo 


’i<ri 


The  cost  of  3  is  defined  to  be  maxlq^}.  We  wish  to  find  a 
schedule  with  minimum  cost. 

As  pointed  out  by  Rinnooy  Kan  [8],  this  problem  may  be 
solved  by  first  determining  tne  finish  time  of  a  minimum 
finish  time  schedule  for  the  n  jobs.  If  the  n  jobs  are 
already  in  nondecreasing  order  of  release  times  r^,  tnen  a 
minimum  finisn  time  schedule  has  start  times  s^  given  by: 

(j.l) 

si»min{ri,  s-^j+l} 


This  equation  is  a  special  case  of  Equation  (6.1)  of 
LiJ  and  so  we  may  use  algorithm  PADMIS  of  Cl]  to  obtain  the 
s-s.  The  finish  time,  C*,  of  a  minimum  finish  time  schedule 
for  the  n  jobs  is  sn  +  1. 

For  simplicity,  let  us  assume  that  r.»  i3.  Once  C*  has 
been  computed,  a  schedule  that  has  minimum  cost  may  be 
obtained  by  constructing  the  convex  bipartite  graph,  <3,  with 
Aa(l,2,...,n}  and  B»{d, 1 , . . . ,C*-1 ) .  There  is  an  edge  from  i 
in  A  to  j  in  B  iff  j>r. .  (A  represents  the  jobs  and  j  *  B 
represents  the  time  slot  Cj/j+1]).  A  cost  c.(j)  is  associ¬ 
ated  with  edge  (i,j). 


From  oar  choice  of  C*,  it  follows  that  the  graph  just 
constructed  has  a  maximum  matching  M  of  size  n.  Every  such 
matching  defines  a  feasible  schedule.  A  matcning  for  which 
the  maximum  edge  cost  is  minimum  defines  the  optimal  match¬ 
ing  for  our  scheduling  problem. 

As  suggested  by  Rinnooy  Kan  [8],  this  matching  can  be 
obtained  by  performing  a  binary  search  on  the  0(nz)  edge 
costs  in  G.  At  each  iteration  of  this  binary  search,  all 
edges  with  cost  greater  than  c  are  deleted  from  G  and  a  max¬ 
imum  matching  in  the  resulting  convex  bipartite  graph  G' 
found.  If  tnis  matcning  is  of  size  less  tnan  n,  then  there 
is  no  matching  in  G  with  maximum  edge  cost  no  more  than  c; 
otherwise  there  is. 


Complexity  Analysis 

o 

The  initial  sort  by  release  times  can  be  done  in  O(log  n) 
time  using  0(n)  PEs  C7].  C*  can  be  found  in  O(logn)  time 
using  n  PEs  and  the  algorithm  of  [lj.  Each  iteration  of  the 
binary  search  requires  the  construction  of  G*  and  the  solu¬ 
tion  of  a2matching  problem.  T£e  first  of  these  can  be  done 
in  J( log  nl  time  using  J(n/logZn)  PEs.  Tne  second  tasx 
taxes  U( log  n)  time  and  0(n)  PEs.  The  total  time  needed  to 
determine  the  matching  with  minimum  cost  is  therefore 
O(log  n)  and  the  number  of  PEs  needed  is  J(n  /log  n) . 
Hence,  the  overall  complexity  of  the  parallel  scneduling 
algorithm  is  G(log  n).  The  number  of  PEs  used  is 
0(n  /log  n) .  The  complexity  of  the  fastest  sequential  algo¬ 
rithm  known  for2this  problem  (i.e.,  the  algorithm  described 
above)  is  G( n^logn)  . Hence,  the  EPU  of  our  parallel  algo¬ 
rithm  is  O(nzlogn/ ( logJn*  nVlog  “)  )=0(  1 )  . 


A’?.  J°k  Sequencing  With  Release  Times  and  Deadlines 

Once  again,  each  job  has  a  processing  requirement  that  is 
one  unit.  However,  this  time  there  is  a  release  time  r., 
deadline  d.,  and  weight  w.  associated  witn  job  i.  Ho  job 
may  be  scheduled  before  its  release  time.  We  are  interested 
in  obtaining  a  maximum  weight  subset  of  the  n  jobs  such  that 
all  jobs  in  this  subset  can  be  scheduled  on  a  single  machine 
in  such  a  manner  that  none  starts  before  its  release  time  or 
finishes  after  its  deadline. 

This  problem  can  be  solved  efficiently  by  a  sequential 
greedy  algorithm.  Jobs  are  considered  in  nondecreasing 
order  of  the  weignts.  If  the  job,  i,  currently  being  con¬ 
sidered  can  be  scheduled  together  with  all  the  others  so  far 
scheduled  such  tnat  no  job  violates  its  release  time  or 
deadline,  then  job  i  is  in  the  optimal  schedule.  Otherwise, 
it  is  not.  This  feasibility  test  can  be  carried  out  in  0(n) 
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time.  So,  the  sequential  algorithm  has  complexity  0(n  ). 

For  the  parallel  algorithm,  we  resort  to  a  convex 
bipartite  graph  formulation.  The  A  3et  represents  the  n 
jobs;  so,  |A|=n.  The  B  set  represents  time  slots.  Let 
u=min(r.j  and  o=max{d.}.  B  =  { a, a+1 , . . . , b-1 } .  In  tne  con¬ 
vex  bipartite  graph  corresponding  to  the  scheduling  problem, 
there  is  an  edge  from  i  in  A  to  j«B  iff  r.£j<d. .  The  weight 
of  each  edge  incident  on  vertex  i*A  is  w. .  tfe^wisn  to  find 
a  matching  d  for  which  the  sum  of  edge  weights  is  maximum. 

our  parallel  algorithm  for  this  is  based  on  the  greedy 
sequential  algorithm  for  joo  sequencing  witn  release  times 
and  deadlines.  Define  the  binary  relation  R  on  tne  weights 
w^  as  follows: 

w^  R  Wj  iff  or  (w^»Wj  and  i>j) 

Let  G^  be  the  convex  bipartite  graph  obtained  from  G  by 
deleting  all  edges  with  weight  w.  such  that  w^  R  w.. 
Clearly,  the  following  determines  whether  or  not  i  is  in  tne 
maximum  weight  matching: 

a)  Determine  the  size  of  the  maximum  cardinality 

matching  in  G. .  Let  this  be  s. 

b)  Delete  all  edges  incident  on  vertex  i-*A  of  G.. 

Determine  the  size  of  the  maximum  cardinality 

matching  in  the  new  graph.  Let  this  be  s'.. 

c)  i  is  in  the  maximum  weight  matching  iff  s^s'^  +  1* 

So,  each  vertex  can  determine,  in  parallel,  whether  or 

not  it  is  in  the  maximum  weight  matching.  To  avoid  read 

conflicts,  we  need  to  raaxe  n  copies  o|  the  graph  G.  This 
can  be  done  in  O(logn)  time  usi^g  0(n  )  PEs.  G^  can  be  con¬ 
structed  in  0(1)  time  using  0(n  ;  PEs  per  vertex  i.  T^en, 
steps  a),  b),  and  c)  above  can  be  performed  in  D(log  n) 

time  using  a  total  of  0(n  /  PEs.  (Note  that  only  pass  1  of 
the  parallel  convex  matching  algorithm  need  be  executed.) 
Tne  total  time  needed  to  determine  the  subset  of  A  in  the 
maximum  weight  matching  is  therefore  O(log  n).  The  number 
of  PEs  can  easily  be  reduced  to  Q(n  /logn) . 

Once  we  have  obtained  the  above  subset  of  A,  the  actual 
matching  can  be  determined  by  deleting  from  G  all  vertices 
that  are  not  in  this  subset.  A  maximum  cardinality  matching 
on  the  resulting  graph  yields  the  desired  matching  with  max¬ 
imum  weight.  This  corresponds  to  an  optimal  schedule  for 
the  scheduling  problem. 

The  complexity  of  the  parallel  algorithm  is  O(log2n) 
and  its  EPU  is  O(nz/(log^n  *  nz/logn))  ■  O(l/logn). 

Exmaple  3^JL:  Figure  3.1  gives  an  example  job  set.  Let 


r* 


-  Id  - 


H.*{jliRjj.  Figure  3.2  gives  the  H.s  for  our  example.  The 
c&nvax  graph  G.  is  shown  in  Figure  3^3.  The  maximum  matcn- 
ing  in  3^  is  3nown  by  broken  lines.  The  maximum  matching  in 
the  graph  obtained  from  3^  by  delating  edges  incident  on 
vertex  4  is  shown  by  wavy^lines.  These  two  matchings  are  of 
the  same  cardinality.  Hence,  vertex  4  is  not  in  the  maximum 
weight  matching.  [] 
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3.3  Preemptive  Scheduling  With  Precedence  Contraints 

bet  J  be  a  3et  of  n  unit  time  jobs.  Let  P  be  a  precedence 
relation  on  J.  For  i«J  and  j*J,  iPj  iff  i  must  be  completed 
before  j  can  commence.  We  may  assume  that  P  is  a  partial 
order.  Hence,  P  may  be  represented  as  a  directed  acyclic 
graph  (dag)  as  in  Figure  3.4  (strictly  speaxing,  P  is  the 
transitive  closure  of  this  graph).  The  directed  edge  <i,j> 
means  that  i  must  be  completed  before  j  can  start  (i.e.,  iPj 
) .  We  also  assume  that  jobs  have  been  indexed  30  that  iP j 
implies  i<j.  This  is  true  of  the  indexing  (of  nodes)  in 
Figure  3.4. 


Muntz  and  Coffman  C 5  3  have  developed  a  level  algoritnm 
to  obtain  minimum  finish  time  preemptive  schedules  for  J  (as 
above)  when  the  number  of  machines  available  is  2.  Their 
algorithm  also  works  for  the  case  when  the  processing  times 
are  mutually  commensurable  (rather  than  simply  1  unit).  A 
set  of  times  (t.,  t^#  ...»  tn)  is  mutually  commensurable  iff 
each  is  a  multiple  of  some"  number  w.  In  the  case  of 
mutually  commensurable  times,  each  job  is  broken  into  a 
chain  of  jobs  each  requiring  w  units  of  processing.  In  this 
section,  we  shall  deal  directly  only  with  the  case  of  unit 
time  tasks. 
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H,  -  {2,3,5,7,3,9,11} 

Si  -  6  Si'-  5  Job  1  is  in 

H2  -  {3,5,7,8,9,11} 

Si  a  5  Si'  -  5  Job  2  is  out 

H3  -  {5,8,11} 

Si  *  4  Si’  -  3  Job  3  Is  in 

H4  -  {1,2,3,5,7,8,9,11} 

Si  -  o  Si'  *  6  Job  4  is  ouc 

Hg  -  {8,11} 

Si  *  3  Si'  -  2  Job  5  is  In 

*  11,2,3,4,5,7,8,9,10,11} 

Si  »  7  Si'  *  6  Job  6  is  In 


H?  -  {3,5,8,11} 

si=5  s,  '*4  Job  7  Is  in 

Hg  a  {11} 

Si  -  2  Si'  -  1  Job  8  is  In 

H9  -  {3,5,7,8,11} 

Si  *  5  Si'  =  5  Job  9  is  out 

H10*  {1,2,3,4,5,7,8,9,11} 

Si  =  6  $i  *  6  Job  10  is  out 

hl1-  {0} 

Si  -  1  Si'  *  0  Job  11  Is  In 


Figure  3.2^ 

Let  3  De  the  dag  representation  of  the  precedence  rela¬ 
tion  P.  A  node  in  G  is  terminal  iff  its  out-degree  is  it). 
The  level  of  a  node  in  G  is  the  length  of  the  longest  path 
from  v  to  any  terminal  node.  All  terminal  nodes  have  level 
•J.  In  Figure  J.4,  the  number  outside  each  node  is  its 
level . 


Muntz  and  Coffman's  £5J  algorithm  to  obtain  minimum 
finish  time  2-machine  schedules  is  given  in  Figure  3.5.  One 
observes  that  the  objective  of  tne  for  loop  of  this  algo¬ 
rithm  is  to  minimize  the  number  of  sets  3.  with  Is. 1*1. 
This  can  be  done  in  parallel  as  follows.  Let  t  be  the  max¬ 
imum  iavsl  in  the  precedence  graph  G.  Let  L^  be  the  level 
of  vertex  i  in  3.  Let  D.  be  the  length  of  theiongest  path 
from  any  vertex  with  no  predecessors  to  i.  Let  H.*L-D. . 
The  CL..D.J  pairs  of  each  vertex  in  our  example  is  given  tn 
Figure  3.4.  tfote  that  H.  is  the  largest  j  such  that  job  i 
could  possibly  be  in  S.  following  the  execution  of  the  for 
loop  in  the  Muntz-Cof fman  algorithm. 

A  job  i  is  said  to  be  critical  iff  L,»H. .  Only  non- 
critical  jobs  are  candidates  £or  displacement1^  other  sets. 
Also,  only  those  sets  S.  that  initially  have  exactly  one 
critical  job  are  candlates  to  receive  a  new  job  during  the 
for  loop.  Uote  that  every  must  have  at  least  one 


2  J 


Figure 

critical  job.  In  Figure  3.6,  the  noncritical  jobs  nave  been 
marked  with  an  X.  Sets  that  are  candidates  to  receive  a  new 
job  have  also  been  marked  with  an  X.  Double  vertical  lines 
mark  sat  boundaries. 

To  determine  the  movement  of  noncritical  jobs  to  candi¬ 
date  sets,  we  construct  a  convex  bipartite  graph  in  which 
the  A  set  consists  of  the  noncritical  jobs  while  the  B  set 
consists  of  the  candidate  sets.  There  is  an  edge  from  i-*A 
to  j«B  iff  L.  <  j  <  ri, .  For  our  example,  the  resulting 
graph  is  shown  in  Fig&re  3.7. 

We  intend  to  use  a  maximum  matching  M  in  the  graph  con¬ 
structed  above  as  follows:  if  (i,j)  is  an  edge  of  the  match¬ 
ing  M,  then  job  i  is  moved  to  set  s.;  vertices  of  A  that  are 
not  matched  to  any  vertex  in  B  remain  in  their  original  set. 

The  maximum  matching  M  is  to  be  obtained  from  the  con¬ 
vex  bipartite  graph  constructed  above  as  follows: 
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Figure  3.4 


procedure  SCH 

l,  <—  maximum  level 

S.  all  jobs  in  J  with  level=i,  £Ki£L 

for  k  «—  L  to  1  do 

if  I 1 =1  then  find  the  largest  j  such  that  j<k  and 
s.  contains  a  job  all  of  whose  pro¬ 
cessors  are  in  S.,  S,  S. 

If  no  such  j  exists  then  exit  enaif . 
Let  q  be  this  value  of  j  and  let  r 
be  the  job  in  S  with  the  above  pro¬ 
perty  ^ 

S  <-  S  -  lr};  S  <-  S  U  {r} 
endif  ^  H 

end for 

preemptively  schedule  the  jobs  in  each  S-  using 
McNaughton's  rule.  1 

Jobs  in  S.  precede  those  in  S.  . . 
end  SCH  1  1-1 


Figure  3^.5^  The  Muntz-Cof fman  Algorithm 
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B 


new  index 


Figure  .3.7,  Convex  Bipartite  Graph. 

(1)  reindex  the  vertices  in  B.  The  new  index  of  vertex 
i«B  is  L-i.  So,  for  i-*A,  s.=L-H.  and  h.=L-L.. 

(2)  use  Glover's  algorithm  considering  vertices  in  B  in 
increasing  order  of  index. 


Let  S'L,S'L  , , ...,S’  be  the  sets  obtained  by  updating 
SL'SL-1' using  the  maximum  matching  above.  These 
updated  sets  satisfy  the  following  properties; 
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(a)  Thera  is  no  joo  pair  (i,j)  such  that  iPj  and  i-*S‘  , 

j-*S'b  for  a  <_b.  a 

(b)  There  is  no  partition  of  the  job  set  J  that  satis¬ 
fies  property  (a)  and  has  a  fewer  number  of  sets  of 

size  1  than  in  S '  r  ,  .  .  . ,  S  ’  ,  . 

L  o 

For  our  example,  the  maximum  matching  M  is 
{ ( 2,4 ) , (4, 2 ) , ( 9 , 1 ) }  (solid  lines  in  Figure  3.7).  So, 
S'.=S.,  i=0, 1,2,5.  S'  ={5,6,7};  and  S'  ={5,2}.  Properties 
(a;  and  (b)  hold  for  our  example. 

The  correctness  of  (a)  follows  from  the  indexing  scneme 
imposed  on  the  jobs  (viz.  iPj  implies  i<j)  and  the  unique 
way  in  which  Glover's  algorithm  works.  Suppose  that  state¬ 
ment  (a)  i3  incorrect.  Then,  there  exists  a  job  pair  (i,j) 
such  that  iPj  and  i*S'a,  j-*S'.D  and  a<b.  Since  no  such  joo 
pair  exists  in  S  ,  ...,  S^,  it  must  be  the  case  tnat  job  j 
is  noncritical  andis  matched  to  b  (old  indexing)  in  M  (as 
only  noncritical  jobs  can  change  sets  and  no  job  moves  to  a 
set  of  smaller  index) .  From  the  definition  of  L  and  H  and 
the  knowledge  that  iPj,  it  follows  that  H.>L^,  d.  _>  L., 
L . >L . ,  and  H->d..  If  job  i  is  critical,  then  it  is  tne  case 
tliat^L  •  >H .  .  1So?  job  j  cannot  possibly  be  matched  to  a  3et  b 
with  aib.  denca,  we  may  assume  that  i  is  also  noncritical. 
Again,  if  i  is  not  matched  in  M,  then  a>b  as  for  a  to  be 
less  than  or  equal  to  b,  it  must  be  the  case  that  L.<H.. 
But  since  h.<h.  and  i  must  get  matched  to  b  befire^j 

can  (see  Glower's  algorithm).  Hence,  i  and  j  must  both  be 
matched  in  M.  i  is  matched  to  a  and  j  to  b.  dence,  a  ^  b. 
If  a<b  then  vertex  b  is  considered  before  a  (as  its  new 
index  is  L-b<L-a).  Since  i  P  b,  it  must  be  the  case  that 
b  >_  L-  •  Also,  L.  >L  ■  implies  that  h^<h..  Hence,  Glover's 
algorithm  will  match  x  to  b  and  not  j  to^b. 

The  truth  of  property  (b)  is  established  using  the  fol¬ 
lowing  reasoning.  If  i*B  is  such  that  no  vertex  of  A  is 
matched  to  i  in  M,  then  Is'. 1=1.  To  see  this,  observe  that 
S.  has  exactly  one  critical  job.  Every  noncritical  job  in 
i3  a  candidate  for  matching  with  i«B.  Thus,  if  no  job  i3 
matched  with  i,  then  all  these  noncritical  jobs  are  matched 
elsewhere  and  I  S'  .  1=1.  It  is  readily  seen  that  if  i<B  is 
matched  in  M,  then  Is'. I  >1.  Hence,  the  number  of  sets  of 
size  1  is  minimized  whenxa  maximum  matching  is  used. 


Complexity  Analysis 

The  3teps  in  our  parallel  preemptive  scheduling  algorithm 
are : 

1)  Determine  L.  and  H.,  l<i£n 

2)  Mark  the  nohcritical  jobs 

3)  Mark  the  candidate  sets 

4)  Construct  the  convex  bipartite  graph 
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5) 

Find  the 

maximum 

3  ) 

Reassign 

matched 

7) 

Scnedule 

jobs  in 

rule . 

matching  M 

joQs  to  chair  new  sets 

the  same  set  using  Mc^a ugh tons 


Step  1  can  be  performed  in  O(log  n)  time  using  the 
critical  path  algorithm  developed  by  Dexel,  tJassimi,  and 
Sanni  L J J •  This  algorithm  uses  0(n  /logn)  PEs .  The  noncrit- 
ical  jobs  can  be  identified  in  0(1)  cime  using  n  PEs.  Tne 
sets  S,,...,S,  may  be  obtained  in  O(logn)  time  by  sorting  on 
L. .  This  requires  0(n  )  PEs.  The  candidate  sets  can  be 
marked  in  O(logn)  time  and  the  resulting  convex  bipartitg 
grapn  may  be  constructed  in  an  additional  0(1)  time  using  ti 
PEs.  The  matching  requires  0(log  n)  time  and  0(n)  PEs  while 
the  reassigning  takes  0(1)  time.  A  parallel  implementation 
of  MCNaughtons  rule  appears  in  LI].  This  has  complexity 
O(logn)  and  uses  0(n/logn)  PEs. 


dance,  the  overall  complexity  of  our  seven  stag  paral¬ 
lel  algorithm  for  preemptive  scheduling  is  O(logn).  The 
number  of  PEs  used  is  0(n/logn).  Since,  the  fastest 
sequential  algorithm  known  for  this  problem  (i.e.,  the 
Muntz-Cof fman  algorithm)  has  complexity  0(n  j,  the  EPU  of 


our  parallel 
0(  1/ ( nlogn) ) . 


algorithm  is 


0(n^/ ( log^n*n^/logn) ) 


4 .  Conclusions 

We  have  developed  fast  parallel  algorithms  for  several  ver¬ 
sions  of  tne  matching  problem  for  convex  bipartite  graons. 
In  Section  2,  we  explicitly  considered  the  maximum  cardinal¬ 
ity  matching  problem.  In  Sections  3.1  and  3.2,  we  con¬ 
sidered  the  problems  of  obtaining  a  maximum  matching  that 
minimized  che  maximum  weight  edge  used  as  well  as  one  tnat 
maximized  the  total  weight  of  the  used  edges.  Both  tnese 
problems  were  discussed  in  connection  witn  the  scheduling 
problems  in  which  they  naturally  arose.  Finally,  the  maximum 
cardinality  matching  algorithm  was  again  used  to  solve  the 
preemptive  scheduling  problem  of  Section  3.3  efficiently. 

This  paper  has  further  enhanced  the  utility  of  the 
binary  tree  metnod  of  Dekel  and  Sahni  [2]  for  the  design  of 
parallel  algorithms.  It  should  also  be  pointed  out  that 
while  all  of  our  complexity  analyses  have  assumed  the  avai¬ 
lability  of  as  many  PEs  as  needed,  our  algorithms  can  be 
used  when  fewer  PEs  are  available.  The  complexity  of  each 
algorithm  will  increase  by  no  more  than  the  shortfall  in 
PEs.  So  if  only  half  the  number  of  PEs  is  available,  then 
the  time  needed  will  at  most  double  (except  for  a  possible 
constant  increase  in  overhead) . 
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